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Abstract 

In the production or annihilation of a pair of fermions, the initial-state or final-state interactions 
often lead to significant effects on the reaction cross sections. For Coulomb-type interactions, the 
Gamow factor has been traditionally used to take into account these effects. However the Gamow 
factor needs to be modified when the magnitude of the coupling constant or the relative velocity 
of two particles increases. We obtain the relativistic generalization of the Gamow factor in terms 
of the overlap of the Feynman amplitude with the relativistic wave function of two fermions with 
an attractive Coulomb-type interaction. An explicit form of the corrective factor is presented for 
the spin-singlet S-wave state. While the corrective factor approaches the Gamow factor in the 
non-relativistic limit, we found that the Gamow factor significantly over-estimates the effects when 
the coupling constant or the velocity is large. 
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I. INTRODUCTION 



The final-state interaction (FSI) and initial-state interaction (ISI) are important processes 
in particle and nuclear physics |l , O, I4 , O, 0, 3, Is > S ll^; III > ll^ • They lead to an 
enhancement of the reaction cross section for attractive interactions and a suppression for 
repulsive interactions. The effects are especially large near the production threshold in the 
region of low-energy annihilation. We shall describe these FSI/ISI in terms of a corrective 
factor which we shall call the ii'-factor. It is defined as the ratio of the cross section with 
the interaction to the corresponding quantity without the interaction. 

In non-relativistic physics the i^-factor can be obtained by solving the two-body 
Schrodinger equation nonperturbatively under their mutual interaction. It is determined 
by calculating the absolute square of the relative wave function at the origin. As is well 
known, for the electric- Coulomb and color-Coulomb interaction, \^(r) = — a/r, the non- 
relativistic corrective i^-factor is the Gamow-Sommerfeld factor P, ^ (or simply called, the 
Gamow factor), 
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where t] is the Sommerfeld parameter. 
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a is the coupling constant (positive for an attractive Coulombic interaction), and v is the 
relative velocity of two particles. For two-equal masses which we shall consider in this paper. 



V is given in terms of their center-of-mass energy ^/s by 
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This gives v ~ 2^/1 — Am'^/s when ^/s ~ 2m and v —>■ 1 when s 00. 

As was pointed out by Chatterjee and her collaborator , the Gamow factor Eq. 

((H), has been traditionally used to study non-relativistic Coulomb-type ISI/FSI in reactions 
involving production or annihilation of particles. The results are then interpolated with the 
well-known perturbative QCD corrective K-factor at high energies, following the procedure 



of Schwinger 



14l |. It predicts an enhancement for qq in color-singlet states and a suppression 



for color-octet states, the effect increasing as the relative velocity decreases. Consequences 



on dilepton production in the quark-gluon plasma, the Drell-Yan process, and heavy quark 
production processes were also examined. 

Although the above Gamow factor gives an approximate description of the FSI/ISI effects 
0], it is useful to obtain a more accurate description as there are physical processes in which 
the interaction coupling constant or the relative velocity of the pair can be quite large and the 
use of the non-relativistic treatment may not be adequate. For example, in the production 
of a charm quark pair, the coupling constant of the color-Coulomb interaction between the 
charm quark and antiquark is about 0.3, which is quite large. Furthermore, as the charm 
quark mass is large, the magnitude of the relative velocity between the produced quark 
and antiquark can be quite small in low-energy production near the threshold. The FSI/ISI 
effects can be quite large for large coupling constants and small relative velocities. Baym and 
P. Braun-Munzinger modified the Gamow factor in their study of the final-state Coulomb 
interaction and effects on the Hanbury-Brown Twiss effects of intensity interferometry . 
A negatively charged particle in a nucleus with a large Z number will also be subject to 
strong FSI/ISI. Although the effect of the interaction is very large for low relative velocities, 
it is nonetheless useful to see how the effect varies as the velocity increases. For brevity 
of notation, we shall use the term "Coulomb interaction" with a variable coupling constant 
to refer to both the electric- Coulomb and color-Coulomb interactions. We shall limit our 
attention to attractive Coulomb interactions, although similar formulation can be carried 
ou. fo. ..epu.s,ve Cou.o.b ,„te..ac«o„. and screened YuUawa „.e.ac«o„s Q. 

Relativistic treatment is needed for strongly attractive interactions, even when the relative 
asymptotic velocity between two particles at r ^ oo is small, as two particles can reach 
relativistic velocities at short distances due to the strongly attractive interaction. Relativistic 
treatment is also needed when the asymptotic relative velocity of two particle approaches 
the speed of light. 

The evaluation of the relativistic corrective if-factor involves the non-perturbative treat- 
ment of the relativistic two-body equation of motion under their mutual interaction. Com- 
pared with the non-relativistic Schrodinger equation involving the Coulomb potential, there 
is an additional attractive effective potential proportional to — |V(r)|^, and a repulsive term 
from the space-like part of the gauge interaction, which lead to a non-trivial behavior when 
the coupling constant becomes large. In the case of fermions with the Coulomb interaction, 
there are further modifications associated with additional spin-dependent potentials. 
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FIG. 1: The Feynman diagrams included in the calculation 



In a previous study [l^ we presented a method to study the relativistic Coulomb FSI/ISI 
effects for a pair of bosons. The corrective i^-factor was evaluated by taking the overlap 
of the relativistic wave function with the corresponding Feynman amplitude. Its analytic 
form was obtained and its numerical values were compared with those of the Gamow factor. 
For attractive interactions, we found that the Gamow factor over-estimates the corrective 
factor for most energies and even more in the relativistic region with a large magnitude of 
the coupling constants. 

We would like to generalize these previous boson results to a pair of fermions with an 
attractive Coulomb-type final state interaction. The fermion results are of more practical 
interest as the electromagnetic gauge field or the chromodynamical gauge field couples to 
fermion fields and the results obtained here may be directly applied to the production or 
annihilation of a pair of fermions in a color-singlet state. A brief report of the present results 
was presented in ref. [l^ 



II. THE K-FACTOR 



We shall be interested in processes involving the production or annihilation of a pair 
of fermions pi and p2 with equal masses by photons (or gluons) of momenta ki and k2 
represented by the Feynman diagrams in Fig. 1. In these diagrams, a solid line represents a 
fermion and a wavy line can be either a photon or a gluon. We shall evaluate the A'-factor 
using a boson- fermion interaction vertex coupling constant g, which will be canceled out in 
the ratio in Eq. (jH} below. The A"-factor depends on the FSI/ISI between the fermions and 
does not depend on how the pair of fermions is produced or annihilated. For definiteness, 
we shall study the production process ki + k2 ^ pi + P2- 

The simplest description of the process is to assume that there is no FSI/ISI and the 
probability amplitude for the production of this pair of particles pi and p2 can be determined 
by means of perturbation theory. The state l^p^pa) P1P2 P^ir after the reaction 



ki + k2 ^ Pi + P2 is represented by the state vector 



l^pipa) = M{kik2 ^PlP2)\piP2), 



(4) 



where A4(fci/c2 P1P2) is the Feynman amphtude for the process. For two-particle system 
P1P2, we define the center-of-mass momentum P = pi + p2 and the relative momentum 
q = {pi — P2)/2. Therefore, pi = P/2 + q and P2 = P/2 — q. 

On the other hand, under their mutual interaction between pi and p2, we can describe 
an interacting piP2 pair with a center-of-mass momentum P as 

|v^y)=V^(q)|P). (5) 

actor defined as the ratio of the cross sections with and 



We introduce the corrective K- 
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where the scalar product (\E'v'|'^'pip2) gives the probabihty amplitude for the produced pair 
I^pjpj) to be in the interacting state l^^v), 



(27r)= 
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and the scalar product (\E'o|'^pip2) is similarly defined in terms of the wave function ipoi'i) 
for a pair of free fermions. The corrective i^T-factor should approach the Gamow factor in 
the non-relativistic limit. 

The cross section with the FSI/ISI corrections is obtained simply by multiplying the 
lowest order cross section with this corrective i^-factor, 



/ Production or annihilation 
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f Production or annihilation 
cross section without FSI/ISI {gq] 



III. DIRAC EQUATION FOR THE COULOMB INTERACTION 



To obtain the two-body wave function , w e use the relativistic two-body equation as 



formulated in Dirac's constraint dynamics 
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23| for two fermions with 



4- momenta pi and p2- We choose to work in the center-of-mass system in which P = (a/s, 0) 



and q = (0, q). In the absence of any interaction, we have the relation between an effective 
energy e^, and a generahzed reduced mass m„, as given by 



0, 



where 



s — 2m^ 



and 



m 
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Next, in the case when the two fermions interact with a mutual Coulomb- type interaction 



V{r) 
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he solution for the two spin-| system under a mutual interaction V{r) can be written as 
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where 



To remove the complications brought by the spinor algebra, we shall carry out calculations 



for the production of the singlet {S = 
following equation of relative motion |l 




svstem. The spin-singlet state is governed by the 



211, 
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By factoring off the angular dependence and the spin dependence: 
Ri{i^)Yimi(^, <P)x''^^^^ ■ The Schrodinger-like radial equation for the state can be written as 
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Ri{r) = 0, 



z^ z z^ 

where z = pr, p is the asymptotic momentum at r ^ oo given by 



(15) 



The wave function Ri{r) can be represented by the dimensionless variable z = pr and is 
characterized by two dimensionless parameters: f] = a/v and a^, where v = p/ew is given 
by Eq. Q. The solution of Eq. (fT3j) is 

Ri{r) = \^e^^/'{2zzr-'/'e-'\F,{a, b, 2iz) (17) 



where 



a = /x+^ + zr/, (18) 



6 = 2^ + 1, (19) 



li=\l{l + \Y-a\ (20) 

and the normalization constant has been determined by using the boundary condition that 
at r ^ oo, Ri{r) — > sin(pr + 5i) /pr with the Coulomb phase shift 5i. For the singlet S-state, 
the critical value of a is 1/2. 

IV. THE FEYNMAN AMPLITUDE AND THE OVERLAP WITH THE 
COULOMB WAVE FUNCTION 

We consider the production of the pair of fermions from the fusion of two gluons (or 
photons) as in the Feynman diagram of Fig. 1. Because the relevant factors associated with 
the mode of production will be canceled out in Eq. © when we take the ratio, the results 
of the ii'-factor depend only on the final-state interaction. 

The diagrams in Fig. 1 give the amplitude 
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where k is the four-momentum of the gluon, q is the four-momentum of one of the fermions, 
and ej is the polarization vector of the i-th gluon. However for the production of a pair 
of fermions under their mutual final-state interactions, we need to project out from the 



Feymann amplitude the proper state iplq) representing the two fermions under final-state 
interactions, 



{-^-M(^i^2 -Pi(q)p2(-q))} - {(2) ^ (1)} 
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Here the factor 1 / and the exchange term is added for the singlet fermion state. Using 
the spinor algebra and full 16 component calculation, the above equation leads to 



where 



2m{E + m) J (27r)3 
J'(q,p, 61,62) 
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(23) 



-i{E + m)^p ■ (61 X 62) - iq ■ 61 p ■ q X 62 
+zq X 61 ■ p q ■ 62 - iq ■ p q ■ 61 X 62. (24) 

As the Coulomb wave function of Eq. ()17j) is given in the configuration space, it is useful 
to write the above integral in terms of the wave function in configuration space. The latter 
is given by 
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In conventional applications, one expands the Eq. 
lowest order g-independent term Aio: 



(25) 

in powers of q and keeps only the 



^A^o + O(|q|). (26) 

In this approximation, Eqs. and give the usual i^-factor as the absolute square of 
the wave function ip{r) at the origin, 

K=\i;{v = 0)\\ (27) 

However, such an approximation cannot be applied to our case with the relativistic wave 
function since the wave function, Eq. (fTTj) . is infinite at the origin. To avoid this singular 
behavior, the full Feynman amplitude is needed to evaluate the overlap integral and the 
i\:-factor in Eqs. (El and ©. 

In terms of the spatial wave function, the overlap integral (|23|1 is 
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We shall specialize to the S-wave case with / = in the present manuscript. Higher partial 
waves can be considered in future work. In this simple case with / = 0, the above integral 
becomes 
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Using the differential equation (fT3j) and the wave function of Eq. (fT7|) . we carry out the 
above integration and obtain 
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6 = m + ip, 
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(32) 



(m + ip)"^ + /c^ 

In Eq. (jHU)) . (a)„ = a(a + l)(a + 2)..(a + n — 1) with (a)o = 1. The quantity (fe)„ is similarly 
defined. This expression looks similar to that in our previous results for boson case in Ref. 



121. 



V. RESULTS FOR THE i^-FACTOR 



We introduce the complex angle variable 

k 



e = tan"^ 
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(33) 



which is a relativistic measure of the relative motion between particles pi and p2. The real 
part of 6 is always tt/A, and the imaginary part is negative, with a magnitude that is half of 
the rapidity of pi(or P2) in the center-of-mass system. 

In terms of the angle variable, Eq.()30|) can be transformed as follows: 
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and V = 3/4 + (/x + ri)/2. To normalize the K-factor, we also need the overlap between the 
Feynman amplitude and the wave function without the final-state interaction. By using the 
wave function ipoi'f') = sin pr/pr for the S-state without the Coulomb potential, we obtained 
the amplitude without the final-state interaction as given by 
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and 0* is a complex conjugate of 0. Then the ratio between the absolute squares of Eqs. ()34|l 
and (jHUj) is the relativistic expression of the /^-factor. 
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We can identify the factor |r(a)e'^''/^/r(6)p as closely related to the Gamow factor G{r]). 
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Therefore, the proper treatment of the dynamics of the interacting particles leads to the 
modification of the Gamow factor G^t]) of Eq. by a factor k, given by 



K = G{ri)n 
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where 
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Here the energy of the fermion is same to the energy of a gluon(i? = fc) in a center of mass 



frame and {E + m)^ + 



2k{k + m). In the limit of a — > or f ^ 0, the factor k 



goes to 1 and the i^-factor is consistent with the Gamow factor. 
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FIG. 2: The ET-factor versus r] for various values of a. 

We note that the center-of-mass energy y/s in units of the rest mass of the produced 
particle is a function of T]/a: 
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(42) 



Various other kinematic variables, such as k/m = ^/s/2m and p/m = A/s74m^^^ can be 
similarly expressed as a function of 7] /a. From these relations and the relation between the 
X-factor and rj and a, we can find out the X-factor for the production of a pair of particles 
in a specific kinematic configuration. 
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FIG. 3: The iC-f actor versus the velocity v for various values of a. 

We showed the behavior of the X-factor as a function of rj in Fig. 2 for various values of 
a. The sohd curve gives the X-factor for a — 0.01 and the dotted curve gives the X-factor 
for a = 0.32. For a fixed value of a, the i^-factor decreases as rj decreases. This is consistent 
with the expectation that the effects of the final-state interaction diminish as the velocity 
becomes relativistic. The limiting value is K—1 as 77 — > Q;(or v ^ 1). Figure 2 also shows 
that for a given value of 77, the X-factor decreases as a increases. It should be noted that 
the same value of 77 corresponds to different velocities v for different values of a. To see the 
effect of final-state interaction as a function of a for a fixed value of we plot in Fig. 3 
the X-factor as a function of v. As one observes, when the velocity is fixed, the X-factor 
increases as the coupling constant increases, indicating a greater effect of the final-state 
interaction as a increases. For all values of a, the i^-factor decreases with v and goes to 
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unity as v approaches 1. The decrease is very rapid for small values of a. 
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FIG. 4: The ratio between the iC-factor and the Gamow factor for various values of a. 

It is of interest to see how the fT-factor obtained here is different from the Gamow factor 
in non-relativistic physics. In Fig. 4(6), we showed the ratio between the fT-factor and the 
Gamow factor for various values of a. As we expect, the ratio is almost 1 for weak coupling 
and the use of the Gamow factor is relatively safe there. However, if we increase a to 0.32, 
the ratio decreases significantly. The Gamow factor over-estimates the magnitude of the 
final-state interaction and therefore it cannot be used for the case with strong coupling. 
There is an effective screening of the long-range Coulomb interaction. As a consequence. 
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the enhancement due to the long-range Coulomb-type interaction is reduced. It can also be 
observed in Fig. 46 that the ratio of K/G{ri) is a relatively slowly varying function of t] for 
1] > 1 but drops down rapidly as 1] decreases in the region of small 1] (high v). 

It is worth pointing out that the expansion of A in Eq. ()35|) is given as a series in powers 
of p/ y6^ + W which increases as the velocity v increases. We still obtain convergent results 
for V up to about 0.8, but there is a limit on using such an expansion for greater velocities 
where p/ is too large to allow for a convergent term-by-term summation. A different 
expansion method may be needed. Fortunately, the i^-factor for this region is so close to 
unity that it can be taken to be unity without incurring much error. 

We can compare the results we have obtained for the case of fermions with those for the 
previous case of bosons. We show in Fig. 4a the ratio of K/G{ri) for the final state interaction 
of two bosons. The K factor for bosons is slightly smaller than the i^-factor for fermions 
and the difference is greater for larger values of a. This means that the overestimations by 
the Gamow factor are larger in bosonic FSI/ISI than in fermionic one. 

VI. CONCLUSIONS AND DISCUSSIONS 

When a pair of particles are subject to final-state interactions, the rate of their production 
is modified. There will be similar effects if the particles interact via initial-state interactions. 
The modification is simplest to be taken into account by using the A'-factor. One calculates 
first the rate for the process when there were no initial- or final-state interactions, using, for 
example, the perturbation theory. The additional initial- or final-state interactions are then 
included by multiplying a A'-factor as given by Eq. (jHl). 

For Coulomb-type interactions, the A'-factor has been traditionally taken to be the 
Gamow factor obtained as the absolute square of the wave function at the origin of the 
relative coordinate. With relativistic Coulomb wave functions, the wave function at the 
origin is infinite and the usual method is not applicable. The A'-factor can be obtained as 
the overlap of the wave function with the Feynman amplitude. 

Our investigation of the A"-factor for the case of the production of a pair of scalar particles 
indicates that there are substantial deviations from the Gamow factor when the strength 
of the coupling is large. In particular, the proper treatment reduces the magnitude of the 
Gamow factor significantly. The reason for this reduction is that in the pair production, 



there is an effective screening of the Coulomb-type interaction arising from the effective 
"exchange" of one of the produced particles. 

We have presented an explicit formula for the relativistic modification of the Gamow 
factor for the production of a pair of fermions. Numerical results are also obtained to show 
the magnitude of the i^'-factor. The results of the i^-factor can be applied to a class of 
processes in which the fermion particles are produced and interacting with a Coulomb- tvpe 



25M 



final-state interaction as for an example in the production of open charm pairs 
Such an application to the production of heavy quarks systems near the threshold will be 
of great interest. 
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